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Abstract 

We discuss various aspects of inclusive top-quark pair production based on Top- 
ixs, a new, flexible program that computes the production cross section at the 
Tevatron and LHC at next-to-next-to-leading logarithmic accuracy in soft and 
Coulomb resummation, including bound-state effects and the complete next-to- 
next-to-leading order result in the qq channel, which has recently become available. 
We present the calculation of the top-pair cross section in pp collisions at 8 TeV 
centre-of-mass energy, as well as the cross sections for hypothetical heavy quarks 
in extensions of the standard model. The dependence on the parton distribution 
input is studied. Further we investigate the impact of LHC top cross section 
measurements at yfs = 7 TeV on global fits of the gluon distribution using the 
NNPDF re-weighting method. 



1 Introduction 



The measurement of the inclusive top-quark pair production cross section is currently of 
great interest at the Large Hadron Collider (LHC). It provides constraints on models that 
explain the anomalous forward-backward asymmetry by new physics; measures the top- 
quark mass cleanly, though not very precisely; is sensitive to the gluon distribution in the 
proton in a momentum-fraction region that is not well constrained by other measurements; 
and it serves as a template for searches for unknown particles with missing transverse 
momentum. For all these reasons a precise cross section calculation is desirable, and, in 
fact, the top-quark pair production cross section is by now the most precisely predicted 
purely hadronic high-energy cross section. 

In [1] some of us presented a calculation of the top-quark pair production cross section 
that included the threshold approximation to the next-to-next-to-leading-order (NNLO) 
partonic cross sections [2], the joint resummation of soft and Coulomb corrections to all 
orders with next-to-next-to-leading logarithmic (NNLL) accuracy, and bound-state effects 
due to Coulomb attraction. The calculation of [lj therefore constitutes the most complete 
implementation of the inclusive production cross section. With the present paper we release 
the first version of the user-friendly public program Topixs ("TOP-pair Inclusive X[cross] 
Section"), which incorporates the calculation of [JJ and an important update as described 
below. The functionality of the program is described in the appendix of this paper. The 
program, together with a detailed manual, can be obtained from the URL 

http : / /users . ph . turn . de/t31sof tware/topixs/ 

which will be continuously updated. Other publicly available implementations beyond NLO 
exist [3H5] and agree at the few-percent level, but do not include soft-gluon resummation 
(HATHOR) and Coulomb resummation (HATHOR, top++) to all orders, while the pro- 
gram provided with j6] obtains the inclusive cross section from resummed or approximated 
differential distributions. Further approximate NNLO results have been obtained in [HE]. 

In the following sections of this work, we consider various topics extending our previous 
results [lj, which are motivated by recent theoretical developments and the anticipated data 
of the 2012 LHC run: 

• We include the (almost) exact NNLO qq partonic cross section, which has recently 
become available |9J. This decreases the cross section in pp collisions at the Tevatron 
by about 1%, together with its theoretical uncertainty, but has a negligible effect on 
pp collisions (LHC). 

• We provide results for yfs = 8 TeV, since meanwhile the LHC centre-of-mass energy 
has been increased to this value. 

• In pQ we exclusively used the MSTW2008 parton distribution functions (PDFs) [TO] . 
Here we compare the predictions for different PDF sets, since, given the accuracy of 
the partonic calculation, the PDF error is now the most important single theoretical 
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uncertainty. Results for different PDFs based on the program HATHOR [3] can also 
be found in [TTj . 

• It is very likely that a perturbative sequential fourth generation will be excluded (or 
discovered) in the 2012 LHC run. Heavy vector-like fermions with Standard Model 
(SM) QCD interactions may exist in other models. We therefore present production 
cross sections for a hypothetical heavy "top" quark. 

• The possibility to constrain the gluon distribution in the proton through the top 
production cross section has been mentioned in [12j[13]. With precise theoretical 
calculations and LHC data at = 7 TeV, this can now be investigated in practice. 
We quantify the impact of the existing data on global fits of the gluon distribution 
using the NNPDF re- weighting method presented in [T^l[T5] . including the case when 
the jet data is left out from the "global" fit. 

The reader interested in the theoretical framework of soft-gluon and Coulomb resummation 
is referred to [T] and jT^ITT]. 



2 Updated cross-sections 



2.1 Theoretical framework for NNLO+NNLL 

The total partonic top-pair production cross section depends on the top-quark mass and 
the kinematic variable /3 = \J\ — Arn^/s. In [1] singular terms as (3 — > were summed to 
all orders at NNLL accuracy by representing the cross section in the form [T7] 



a 



NNLL/ 
pp' \ 



H*,(m t ,ii) I du J R (E - |) W R ( 



(2.1) 



R=l,8 



where E = y§ — 2m t is the energy relative to the production threshold. The soft func- 
tions Wr sum logarithms of /3 from soft-gluon emission. The potential function Jr sums 
Coulomb-gluon exchange related to the attractive or repulsive Coulomb force in the colour- 
singlet and octet channels, respectively. Logarithms of /3 from non-relativistic physics are 
currently included at NNLO in Jr. The hard functions Hr account for the short-distance 
production process of the tt pair. The momentum-space approach to threshold resumma- 
tion [T8l[T9] is employed, where renormalization-group equations evolve the hard functions 
from a hard scale \Xh to the factorization scale \if and the soft function from the soft scale 
H & ~ m t (3 2 to p,f. As default the values fih = 2m t , /// = m t and a "running" soft scale 
Us = 2m t max[/9 2 , f3^ ut ] is used. An alternative method with a fixed soft scale, and the 
determination of (3 CVLt are discussed in detail in pp. 

Since (12. ip is strictly valid only for small /3, the expression is matched to the available 
fixed-order computations. In [TJ the final results were based on the matching equation 
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where 2 \s) is (12. ip expanded to 0{a 2 s ) relative to the Born cross section, and the last 

two terms represent the fixed-order NLO cross section plus the threshold approximation of 
the NNLO terms [[2J. That is, the partonic cross section consisted of all known fixed-order 
terms with the 0(a™) terms with n > 2 from the resummation formula (12. ip addedJl] 

This can now be improved. The full NNLO partonic cross section, not restricted to 
small /3, has recently been computed [9 J for the qq initial state, except for the partonic 
processes qq — > tiqq with same-flavour (anti-)quarks in the initial and final state, which 
are presumed to be small. The matching equation (12.21) is straightforwardly modified to 
include this important result by substituting 

^(*) + *S££(S) ST°(*) (2-3) 

for pp' = qq, where <3"^y fLO (s) is simply the full fixed-order NNLO cross section as given 
in [9], now including all lower orders, and the full renormalization and factorization scale 
dependence [20] @ In the following this implementation will be used and denoted NNLL. 
Note, however, that we keep the labels NNLO app and NNLL 2 for the LHC even though 
the full qq channel is included as well, since at the LHC the gg channel dominates and the 
effect of the replacement in (12.31) is very small. 

The theoretical uncertainties will be computed as detailed in [1] with two modifications: 
(1) In the qq channel the theoretical uncertainty associated with the unknown NNLO result 
parameterized by the estimate of a constant term in the threshold expansion has become 
irrelevant. This uncertainty is now included only for the gg channel. (2) Since we now 
include different PDF sets, we present the PDF+a s error at the 68% rather than 90% 
confidence level (CL). 



2.2 Tevatron results 

Since top pairs in pp collisions at the Tevatron are produced predominantly through the 
qq initial state, we are now in the position to present a (nearly complete) NNLL result 
matched to the (nearly complete) NNLO fixed-order calculation. We update our results 
from p] using the same top-quark pole mass m ( = 173.3 GeV and MSTW2008NNLO [10] 
PDF input with a s (Mz) = 0.117lio!ooi4 EI]- We compare the various approximations in 
Table [H The NLO, NNLO app and NNLL2 entries are the same as in [1] except for the 
PDF+a s error which is now given at 68% rather than 90% CL. With the new NNLO 
result for the qq channel included (second line in the Table), the cross section is reduced 
by about l%o The theoretical uncertainty of the resummed cross section from scale and 

1 This result was denoted NNLL2 in pp. A second matching option, denoted NNLLi, was also considered. 
Since this implementation is no longer useful, once the exact NNLO result is known, we do not discuss it 
in the present paper. 

2 Note that in the expressions for and /|| in Eqs. (Al) and (A2) of [50], the functions ff^ multiplying 
the n ^ 1 terms must be replaced by their threshold limit. 

3 Note that the difference between NNLO and NNLO app on the one hand, and NNLL and NNLL 2 on 
the other is not exactly the same, since we recompute the value of /9 cu t in the resummation method 2, 
see [I], which changes from 0.35 for NNLL 2 to 0.38 for NNLL. 
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NLO 


NNLO a pp 


NNLL 2 


a co +0.36+0.23 
U.UO _o.75-0.22 


7 nfi +0.24+0.10+0.29 
' - uu -0.33-0.10-0.24 


7 99+0.29+0.10+0.30 
' -0.46-0.10-0.25 




7 nn +0.21+0.02+0.29 
' - uu -0.31-0.02-0.25 


7 i c +0.21+0.02+0.30 
' - l ° -0.20-0.02-0.25 




NNLO 


NNLL 



Table 1: The total top-pair cross section (in pb) for m t = 173.3 GeV in pp collisions at the 
Tevatron (y/s = 1.96 TeV) with MSTW2008NNLO (NLO for NLO) input. The first set of errors 
refers to scale variation (scale variation+resummation ambiguities for NNLL 2 and NNLL), the 
last to the 68% CL PDF+o; s error. The second set of errors for NNLO( app )/NNLL(2) arises from 
variations of the unknown NNLO constant term. 



resummation ambiguities and the unknown constant is reduced from ±5.4% to ±2.9%. The 
reduction of the error from NNLO app to NNLO is smaller which we attribute to the fact 
that a non-negligible part of the scale uncertainty originates from the gluon-gluon initial 
state (see Table 5 in p]) so a further reduction is expected from the full NNLO calculation 
for all initial states. The smallness of the correction shows a posteriori that the threshold 
approximation works rather well (though, partly, by a fortuitous cancellation in the large 
(3 region, see [9]). Our best result for the Tevatron tt production cross section, assuming 
m t = 173.3 GeV and a s (M z ) = 0.1171, is therefore 

aa = 7.15 ±g:|J (theory) +° f 5 (PDF±a s ) pb (MSTW2008NNLO). (2.4) 

The cross section increases to 7.26 pb, see Table [6] below, when the "world average" value 
a s (Mz) = 0.1180 is used. We estimate the theory uncertainty by adding the errors from 
scale variation, resummation ambiguities and the estimate of the constant NNLO term in 
the gluon channel in quadrature, where the latter has negligible effect at the Tevatron. 
Table [2] provides the corresponding results in the range mt = 165-180 GeV. 

It is instructive to compare the effect of adding the full NNLO qq partonic cross section 
to the previous NNLL resummations in the momentum-space or SCET-based resummation 
formalism employed here and in p], and to the Mellin-space formalism [HE]. This is 
particularly so as the Mellin-space result as well as the result from [6J have been somewhat 
lower than those given in [TJ, see the comparison [22J. For the central values, we find 

7.22pb([lJ) -»■ 7.15pb (Eq. O) (momentum-space), 
6.72pb(01) -»■ 7.07 pb ([9]) (Mellin-space), 

for identical m t , a s and PDF input. We observe that the two approaches are now in 
full agreement, and that the momentum-space result has turned out to be substantially 
less sensitive to the inclusion of the full NNLO qq channel than the Mellin-space ond_|. 

4 This is partly due to the change in the matching procedure in the Mellin-space results (from 
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Table 2: Total cross sections in pb at the Tevatron for m t = 165-180 GeV. See caption of 
Tabled] for the definition of the various errors. 



m t [GeV] NLO 



NNLO 



NNLL 



1 p c 

loo 


8.70 


+0.48+0.30 
-0.98-0.30 


9.10 


+0.27+0.04+0.38 
-0.40-0.04-0.32 


9.29 


+0.29+0.04+0.40 
-0.26-0.04-0.33 


1 PP 

loo 


O A O 

8.42 


+0.46+0.29 
-0.95-0.29 


8.82 


+0.26+0.03+0.37 
-0.38-0.03-0.31 


n c\c\ 

9.00 


+0.28+0.03+0.38 
-0.25-0.03-0.32 




8.1b 


+0.45+0.28 
-0.92-0.28 


O C A 

8.54 


+0.25+0.03+0.36 
-0.37-0.03-0.30 


8.71 


+0.27+0.03+0.37 
-0.24-0.03-0.31 


loo 


( .y(J 


+0.43+0.27 
-0.89-0.27 


o.Zl 


+0.25+0.03+0.35 
-0.36-0.03-0.29 


Q A A 

8.44 


+0.26+0.03+0.36 
-0.24-0.03-0.30 


169 


7.65 


+0.42+0.26 
-0.86-0.26 


8.01 


+0.24+0.03+0.33 
-0.35-0.03-0.28 


8.18 


+0.25+0.03+0.34 
-0.23-0.03-0.29 


170 


7.41 


+0.40+0.26 
-0.84-0.25 


1 .lb 


+0.23+0.03+0.32 
-0.34-0.03-0.27 


7.92 


+0.24+0.03+0.33 
-0.22-0.03-0.28 


171 


7.18 


+0.39+0.25 
-0.81-0.24 


7.52 


+0.22+0.03+0.31 
-0.33-0.03-0.26 


'IPO 

7.68 


+0.23+0.03+0.32 
-0.22-0.03-0.27 




P. QtP. 

o.yo 


+0.38+0.24 
-0.78-0.23 


7 9Q 


+0.22+0.03+0.30 
-0.32-0.03-0.25 


7 /l/l 


+0.22+0.03+0.31 
-0.21-0.03-0.26 


173 


6.74 


+0.37+0.23 
-0.76-0.23 


7.07 


+0.21+0.02+0.29 
-0.31-0.02-0.24 


7.21 


+0.21+0.02+0.30 
-0.20-0.02-0.25 


174 


6.54 


+0.35+0.23 
-0.74-0.22 


6.85 


+0.20+0.02+0.28 
-0.30-0.02-0.23 


6.99 


+0.21+0.02+0.29 
-0.20-0.02-0.24 


175 


6.34 


+0.34+0.22 
-0.71-0.21 


6.64 


+0.20+0.02+0.27 
-0.29-0.02-0.22 


6.78 


+0.20+0.02+0.28 
-0.19-0.02-0.23 


176 


6.14 


+0.33+0.21 
-0.69-0.20 


6.44 


+0.19+0.02+0.26 
-0.29-0.02-0.22 


6.57 


+0.19+0.02+0.27 
-0.19-0.02-0.22 


177 


5.96 


+0.32+0.21 
-0.67-0.20 


6.24 


+0.19+0.02+0.25 
-0.28-0.02-0.21 


6.38 


+0.19+0.02+0.26 
-0.18-0.02-0.22 


178 


5.78 


+0.31+0.20 
-0.65-0.19 


6.06 


+0.18+0.02+0.25 
-0.27-0.02-0.20 


6.18 


+0.18+0.02+0.25 
-0.18-0.02-0.21 


179 


5.60 


+0.30+0.20 
-0.63-0.18 


5.87 


+0.17+0.02+0.24 
-0.26-0.02-0.20 


6.00 


+0.17+0.02+0.25 
-0.17-0.02-0.20 


180 


5.43 


+0.29+0.19 
-0.61-0.18 


5.70 


+0.17+0.02+0.23 
-0.25-0.02-0.19 


5.82 


+0.17+0.02+0.24 
-0.17-0.02-0.20 



Furthermore, the comparison to the full NNLO result shows that approximate NNLO 
expressions obtained from integrating approximate invariant mass or one-particle inclusive 
distributions |6] do not lead to a better approximation than the threshold expansion. 

It is interesting to study the impact of the new, almost complete, NNLO+NNLL pre- 
dictions on the extraction of the pole mass of the top quark from the measured total 
cross section at the Tevatron. The DO collaboration [23] extracted the pole and MS mass 
using several higher-order predictions of the cross section. For instance, the pole mass 
m t = 167.5^4 7 GeV was obtained using the approximate NNLO calculation from [20] and 
MSTW2008 PDFs. We now study the impact of our updated theory predictions on the 
mass measurement. To this end we fit the NNLL cross section in Table [2] by a function of 

NLO+NNLL in 4 to NNLO+NNLL in 19 ) and the different treatment of the NNLO constant term 
in [TJ and 0]. 
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the form 



/172 5\ 4 

0tfK)= f J (co + CxK- 172.5) + c 2 K-172.5) 2 + c 3 (m i -172.5) 3 ) pb , 

(2.5) 

where all the masses are given in GeV. For the coefficients in the fit ( 12. 5ft we find Co = 
7.325 ± 0.219 ± 0.304, Cl = -(5.8651 ± 0.3163 ± 0.3382) x 10~ 2 , c 2 = (2.4884 ± 0.7131 ± 
0.9027) x 10" 4 , c 3 = (3.7436 ± 4.1254 ± 1.9924) x 10~ 6 where the first error denotes 
the total theory error and the second the PDF+a s error. The maximum of the positive 
and negative errors has been used in the fit. An analogous fit for the dependence of the 
experimental cross section on the reference mass is given by the DO collaboration in [21] for 
the combination of the dilepton and lepton+jet results that corresponds to a t t — 7.56loi56 
pb for a reference mass of wi[ = 172.5 GeV. As in [1] the most probable value of the 
pole mass is obtained by the maximization of a joined likelihood function obtained from 
a convolution of two normalized Gaussians centered at the theoretical prediction and the 
experimental mean value, respectively, and with a width given by the corresponding errors. 
This method is similar to the one used in [23], however, we do not introduce separate 
likelihood functions for scale and PDF uncertainty, but add theory and PDF errors linearly. 
Using our best NNLL prediction, we obtain^ 

m t = 171.4lf;} GeV. (2.6) 

This result is in good agreement with the combined result m t = 173.2±0.8 GeV from the di- 
rect mass determination at the Tevatron [25] as well as our previous result [1] m t = 169.8l 4 'y 
GeV using an ATLAS cross-section measurement [26J and the value m t = 170.3lg'y GeV 
obtained from a CMS cross-section measurement [27J using the calculation of [20J . In ( 12. 6 j) 
the mass used in the kinematic reconstruction of the top-events is identified with the pole 
mass. Allowing for a difference of ±1 GeV in the two mass definitions leads to an addi- 
tional uncertainty of 0.5 GeV. To study the effect of the exact NNLO result of [9] and 
of NNLL resummation, we have repeated the same procedure for the approximate NNLO 
results presented in [T] and the NNLO results in Table HJ We obtain m t = 171.01U for the 
approximate NNLO calculation and m t = 170.5lg4 f° r ^he NNLO result, respectively. 
Consistent with Tabled! the effect of the result of [9] on the NNLO calculation is moderate, 
while a stronger reduction of the error is observed for the NNLL prediction. 



2.3 LHC results 

In this section we provide predictions for the top-antitop total cross section, computed 
with Topixs, for the LHC with the current centre-of-mass energy of 8 TeV, extending 

J The uncertainty is not improved compared to :23] that used selection criteria resulting in a weaker 
TOt-dependence of the experimental cross section than in [24], but a larger cross section measurement of 
<j t t = 8.13lJ'go pb for a reference mass of mt — 172.5 GeV. Naively rescaling the mass-dependence of 
the experimental cross section in [23] to the combined central value and error of [24j we obtain m t = 
171. 6l 4 3 GeV which might indicate the precision that could be achieved with the Tevatron data. 
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Table 3: Total cross sections in pb at the LHC (y/s = 7 TeV) for m t = 165-180 GeV. See 
caption of Tabled] for the definition of the various errors. 



-23.4-6.8 



m t [GeV] NLO NNLO app NNLL 2 

+25.5+8.6 9n7 a +13.6+5.7+9.2 9flQ a +6.6+5.7+9.3 

-27.4-7.8 zu ' -^-14.0-5.7-8.5 ZUJ -^ -6.9-5.7-8.6 

+24.7+8.4 9m -,+13.1+5.5+8.9 9n9 Q +6.4+5.5+9.0 

-26.5-7.6 ^Ul.l -13.6-5.5-8.2 ZUZ - J -6.7-5.5-8.3 

+23.9+8.1 i Qf- n +12.7+5.3+8.6 i Qfi o +6.3+5.3+8.7 

-25.7-7.4 1JO.U -13.1-5.3-8.0 iJU -° -6.5-5.3-8.1 

+23.1+7.9 -i on n +12.2+5.1+8.4 i Qn o +6.1+5.1+8.5 

-24.9-7.2 lOJ.U -12.7-5.1-7.8 lyu -° -6.3-5.1-7.9 

+22.4+7.7 i oq q +11.8+5.0+8.1 i oc n +5.9+5.0+8.3 

-24.2-7.0 ioo -° -12.3-5.0-7.5 lOJ.U -6.1-5.0-7.6 

1 77 8 +H-4+4.8+7.9 -i 7Q r +5.7+4.8+8.0 

1 ' ' -° -11.9-4.8-7.3 1 ' y ' d -5.9-4.8-7.4 

+21.0+7.2 -I 79 r +11.0+4.6+7.7 -i 7/1 i +5.5+4.6+7.8 

-22.7-6.6 1 -11.5-4.6-7.1 1 ' ^- 1 -5.8-4.6-7.2 

+20.3+7.0 -i fi7 /, +10.7+4.5+7.5 -i fi o Q +5.4+4.5+7.6 

-22.0-6.4 iu ' -^-11.1-4.5-6.9 iuo - y -5.6-4.5-7.0 

+ 19.7+6.8 -1 fi9 j +10.3+4.3+7.3 -i fi o q +5.2+4.3+7.4 

-21.4-6.2 i-VZ.t -10.8-4.3-6.7 iuo - y -5.4-4.3-6.8 

+ 19.1+6.6 i p- 7 fi +10.0+4.2+7.1 i cq 1 +5.1+4.2+7.3 

-20.7-6.1 LO ' - u -10.5-4.2-6.5 ±,JJ - 1 -5.3-4.2-6.6 

+ 18.5+6.5 i e-o n +9.7+4.1+6.9 i fa a +4.9+4.1+7.0 

-20.1-5.9 lOO.U -10.1-4.1-6.4 i,J ^-^ -5.1-4.1-6.4 

+ 17.9+6.3 i « a r +9.4+3.9+6.7 1/IQ n +5.3+3.9+6.8 

-19.5-5.7 l^o.O -g.8-3.9-6.2 iq:J - J -5.0-3.9-6.3 

+ 17.4+6.1 1/1/10 +9-1+3.8+6.5 i j r r +4.6+3.8+6.6 

-19.0-5.6 1^.^ -9.5-3 8-6.0 i ^' J -' J -4.9-3.8-6.1 

+ 16.9+5.9 i /in n +8.8+3.7+6.3 i /.i o +4.5+3.7+6.4 

-18.4-5.4 l^U.U _g 2_3.7_5.9 J-^ 1 - -4.7-3.7-5.9 

+ 16.3+5.8 i qc n +8.5+3.6+6.1 i q 7 q +4.4+3.6+6.2 

-17.9-5.3 lOU.U -8.9-3.6-5.7 io ' ° -4.6-3.6-5.8 

+ 15.9+5.6 i q 9 i +8.2+3.5+6.0 i qq q +4.3+3.5+6.1 

-17.3-5.1 IOZ.1 _ 8 . 7 _3 i5 _ 5 . 5 J- 00 - -4.5-3.5-5.6 



165 


203.9 


166 


197.7 


167 


191.6 


168 


185.8 


169 


180.2 


170 


174.7 


171 


169.5 


172 


164.4 


173 


159.6 


174 


154.8 


175 


150.3 


176 


145.9 


177 


141.6 


178 


137.5 


179 


133.5 


180 


129.7 



the results for 7 and 14 TeV given in [T] (the result for m t = 173.3 GeV with the setup 
of [1] was already presented in [28J). We also update the results for 7 TeV centre-of-mass 
energy by including the full NNLO result for the gg-initial state, although the effect on the 
central value is negligible since ti production at the LHC is dominated by the gluon-fusion 
channel where the full NNLO corrections are not known yet. Since the LHC centre-of-mass 
energy after the 2013-2014 shutdown has not been fixed, we do not give results for higher 
centre-of-mass energies. Results for other centre-of-mass energies can be easily obtained 
with Topixs since \fs can be freely chosen in the configuration file, see the Appendix. 
Our best result for the ti production cross section at the LHC, using m t = 173.3 GeV, the 
MSTW2008 PDF sets [TO] and a s {M z ) = 0.1171, is given by 

u t i (LHC, 7 TeV) = 162.4 (theory) _£f (PDF+a s ) pb (MSTW2008NNLO) (2.7) 

<x tf -(LHC, 8 TeV) = 231.8 (theory) (PDF+« S ) pb (MSTW2008NNLO) (2.8) 
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Table 4: Total cross sections in pb at the LHC (y/s = 8 TeV) for m t = 165-180 GeV. See 
caption of Tabled] for the definition of the various errors. 



mt [GeV] 


NLO 




NNLO 


app 


NNLL 2 


165 


290.1 


+36.1+11.7 
-38.1-10.4 


294.5 


+20.0+8.3+12.3 
-20.1-8.3-11.4 


9Q7 n+ 9 - 1 + 8 - 3 + 12 - 4 

Z» 1 .U_g g_g 3_u 5 


166 


281.4 


+35.0+11.4 
-37.0-10.1 


285.7 


+19.4+8.0+12.0 
-19.5-8.0-11.0 


988 1 +8.9+8.0+12.1 
ZOO.±_ 9 3 _ 8 _ n 2 


167 


273.0 


+33.9+11.1 
-35.9-9.8 


277.2 


+18.7+7.8+11.6 
-18.8-7.8-10.7 


97Q c+8. 6+7.8+11. 8 
z ' J -°-9. 1-7.8-10. 9 


168 


264.9 


+32.8+10.8 
-34.8-9.6 


269.0 


+18.1+7.5+11.2 
-18.3-7.5-10.4 


971 9+8.4+7.5+11.4 
z ' i - z -8. 8-7.5-10. 6 


169 


257.0 


+31.8+10.4 
-33.7-9.3 


261.0 


+17.5+7.3+11.0 
-17.7-7.3-10.2 


ocq 9+8.1+7.3+11.1 
zuo - z -8. 6-7.3-10. 3 


170 


249.5 


+30.8+10.2 
-32.7-9.0 


253.3 


+17.0+7.0+10.7 
-17.1-7.0-9.9 


9 p;p; cr+7.9+7.0+10.8 
zdd - d -8. 3-7.0-10.0 


171 


242.1 


+29.9+9.9 
-31.8-8.8 


245.9 


+16.4+6.8+10.4 
-16.6-6.8-9.6 


9/18 n+ 7 - 6 + 6 - 8 + 10 - 5 

z ^°- u -8. 1-6. 8-9.7 


172 


235.1 


+29.0+9.6 
-30.8-8.6 


238.8 


+15.9+6.6+10.1 
-16.1-6.6-9.3 


9AD O+7.4+6.6+10.2 
z ^ u -°-7.9-6. 6-9.4 


173 


228.3 


+28.1+9.3 
-29.9-8.3 


231.9 


+15.4+6.4+9.8 
-15.6-6.4-9.1 


ooq q+7.2+6. 4+9.9 
^°°- y -7.6-6. 4-9.2 


174 


221.7 


+27.3+9.1 
-29.1-8.1 


225.2 


+14.9+6.2+9.5 
-15.1-6.2-8.8 


997 1 +7.0+6.2+9.7 
zz ' • ± -7.4-6. 2-8.9 


175 


215.3 


+26.4+8.8 
-28.2-7.9 


218.7 


+14.4+6.0+9.3 
-14.6-6.0-8.6 


99D «+6. 8+6. 0+9.4 
zzu.U_ 72 _ 6 .0-8.7 


176 


209.1 


+25.6+8.6 
-27.4-7.7 


212.5 


+14.0+5.8+9.0 
-14.2-5.8-8.4 


91/1 q+6. 6+5. 8+9.1 
z± ^-°-7.0-5. 8-8.5 


177 


203.2 


+24.9+8.3 
-26.6-7.5 


206.4 


+13.5+5.6+8.8 
-13.8-5.6-8.1 


908 9+6.4+5.6+8.9 
zuo - z -6. 8-5. 6-8.2 


178 


197.4 


+24.1+8.1 
-25.9-7.3 


200.6 


+13.1+5.5+8.6 
-13.3-5.5-7.9 


909 q+6. 2+5. 5+8.7 
^ uz -°-6. 7-5. 5-8.0 


179 


191.8 


+23.4+7.9 
-25.1-7.1 


194.9 


+12.7+5.3+8.3 
-12.9-5.3-7.7 


196 R+6-1+5.3+8.4 
lyu - u -6. 5-5. 3-7.8 


180 


186.4 


+22.7+7.7 
-24.4-6.9 


189.5 


+12.3+5.1+8.1 
-12.6-5.1-7.5 


1Q1 1+5.9+5.1+8.2 
Ayl -- L -6.3-5.1-7.6 



The NLO, NNLO and NNLL tt cross sections at the LHC with y/s = 7, 8 TeV for the mass 
range m t = 165-180 GeV are shown in Table[3]and|H respectively. The effect of corrections 
beyond the next-to-leading order at the LHC with y/s = 8 TeV is similar to that observed 
for \fs = 7 TeV, increasing the cross section by 2.4% compared to the NLO result, with 
fixed-order NNLO contributions accounting for 1.6%. It can be seen from Tables [3] and H] 
that the uncertainty from the unknown fixed-order NNLO corrections for the gg (and qg) 
initial states, as estimated by the variation of the constant term in the threshold expansion, 
are of the same order as the remaining theory errors of the NNLL results at the LHC. This 
source of uncertainty will be removed once the full NNLO calculation has been performed 
also for these channel^. 

Comparing to our previous predictions for the LHC with y/s = 7 TeV [lj and y/s = 

6 For recent work constraining the NNLO cross section using information on the (3 — > 1 or large invariant- 
mass limits, see [291130]. 
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Table 5: Values for /3 cut . The strong coupling constant is a. s [Mz) = 0.118. 







Tevatron 






LHC7 






LHC8 




PDF set 


k s — 


1 2 


4 


k s = 1 


2 


4 


k s = 1 


2 


4 


MSTW2008 


0.56 


0.38 


0.34 


0.74 


0.54 


0.33 


0.76 


0.54 


0.33 


NNPDF2.1 


0.53 


0.29 


0.19 


0.74 


0.51 


0.25 


0.75 


0.52 


0.25 


ABM11 


0.54 


0.36 


0.22 


0.73 


0.54 


0.33 


0.74 


0.54 


0.33 


CT10 


0.57 


0.39 


0.24 


0.74 


0.54 


0.33 


0.75 


0.54 


0.33 



8 TeV [2S] it is seen that in our approach the inclusion of the exact NNLO result for the 
qq channel amounts to a minimal change of the central value by —(0.1-0.2) pb: 

162.6 pb (|H) -> 162.4 pb (Eq. (E2D) (7 TeV), 
231.9pb ([28]) -> 231.8 pb (Eq. (E3D) (8 TeV). 

The effect on the theory error is somewhat larger and we find a reduction from l^g to 
t 6 6 7 9 pb at y/s = 7 TeV and from t\oi to -li pb at ^/i = 8 TeV. Again, a larger effect on 
the central value is observed in the approach of [U|3]0 

158.7pb([4j) 156.2pb(|5]) (7 TeV), 
226.6pb([4j) 222.5pb(|5]) (8 TeV). 

While these results agree with ours within the quoted uncertainties, the relative difference 
is larger than for the Tevatron which indicates the relevance of the uncalculated fixed-order 
NNLO terms for the dominant gluon-fusion channel. 

3 Dependence on PDF sets 

In this section we present results for the cross section computed with different PDF sets. As 
representatives we choose MSTW2008 [10], NNPDF2.1 [31], ABM11 rjJJ, and CT10 [32]. 
(Note that CT10 at the moment provides only NLO PDFs.) We use m t = 173.3 GeV, and 
choose PDF sets for the common value a s (Mz) = 0.118 of the strong coupling constant. 
This provides complementary information to [llj, where results for different PDF sets 
with their best, but often rather different a s values have been presented, which, due to the 
dependence of the partonic cross sections on a s , obscures the difference due to the parton 
densities themselves. 

7 When comparing the NNLO app results from Topixs to those obtained with top++1.2 [5] one finds a 
small difference of 0.1-0.3 pb. This is due to the implementation of the approximate NNLO corrections for 
the gg initial state, where the exact tree cross sections for the separate colour-singlet and octet channels 
are used in Topixs, while the colour averaged tree cross sections are used in top++. In the threshold 
limit both implementations agree, so they are formally of the same accuracy. 
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Table 6: Total cross sections in pb at the Tevatron and LHC (yfs = 7TeV and 8TeV for 
different PDF sets and common strong coupling input a s (M z ) = 0.118. See caption of Tabled] 
for the definition of the various errors. 



Tevatron NLO NNLO NNLL 



MSTW2008 


6.48 


+0.334 
-0.70- 


-0.23 
-0.22 


7.11 


+0.22+0.03+0.29 
-0.32-0.03-0.24 


7 9 £ +0.22+0.03+0.30 
' - zu -0.21-0.03-0.25 


NNPDF2.1 


6.77 


+0.344 
-0.62- 


-0.23 
-0.23 


7.18 


+0.22+0.02+0.26 
-0.17-0.02-0.26 


7 AQ +0.26+0.02+0.27 
' -0.25-0.02-0.27 


ABM11 


6.45 


+0.264 
-0.64- 


-0.14 
-0.14 


7.29 


+0.22+0.01+0.15 
-0.31-0.01-0.15 


7 a c +0.22+0.01+0.16 
' -0.20-0.01-0.16 


CT10 


6.75 


+0.354 
-0.74- 


-0.46 
-0.35 


7.33 


+0.22+0.03+0.52 
-0.33-0.03-0.40 


7 47+0.23+0.03+0.52 
' -0.22-0.03-0.41 



LHC7 NLO NNLO a pp NNLL 2 



MSTW2008 


151.8 


+18.2+6.8 
-19.9-6.2 


164.3 


+ 10.6 
-11.1 


+4.4+7.2 
-4.4-6.7 


165.9 


+5.3+4.4+7.3 
-5.6-4.4-6.8 


NNPDF2.1 


157.0 


+19.1+6.8 
-19.7-6.8 


161.9 


+8.3+4.3+7.1 
-7.5-4.3-7.1 


166.1 


+7.0+4.3+7.3 
-8.3-4.3-7.3 


ABM11 


123.3 


+14.5+5.0 
-16.1-5.0 


146.7 


+9.4+3.8+5.8 
-9.7-3.8-5.8 


148.2 


+5.3+3.8+5.9 
-5.3-3.8-5.9 


CT10 


148.7 


+17.7+12.7 
-19.3-11.4 


160.7 


+ 10.4+4.4+13.6 
-10.8-4.4-12.2 


162.3 


+5.2+4.4+13.7 
-5.4-4.4-12.3 



LHC8 



NLO 



NNLO 



app 



NNLL 2 



MSTW2008 
NNPDF2.1 
ABM11 
CT10 



91 7 c- +26.0+9.2 
' >J _27.9-8.3 

oop; -i +27.3+9.0 
Jl —27.6—9.0 

1«n o +21.4+6.8 
18U.8_ 23 3 _ 68 

oi 9 o +25.3+16.5 
-27.1-14.8 



no a a +15.7+6.5+9.7 
zo ^-^ -16.0-6.5-9.0 

909 9+12.2+6.5+9.4 
ZOZ.Z -n.o-6.5-9.4 

91 9 1 +14.2+5.7+7.8 
^-L^-J- -14.3-5.7-7.8 

99Q 9+15.5+6.4+17.6 
-15.6-6.4-15.8 



236.5 



+7.4+6.5+9.8 
-7.3-6.5-9.1 



900 n +9.7+6.5+9.7 

ZOO.U -n.7-6.5-9.7 

91 a 1 +7.4+5.7+7.9 
^-L^- 1 -7.2-5.7-7.9 

901 9+7.2+6.4+17.7 
zo± - z -7.2-6.4-15.9 



For NNPDF2.1 the total number of replicas used is A^ rep = 100. To determine the 
PDF+a: s error, this number is distributed over sets with different a s values according to 
a Gaussian distribution with a standard deviation of 5a s = 0.0015. For MSTW2008 and 
CT10 we provide asymmetric PDF+a s errors using the methods from [2"TII3"3"] . In the latter 
case we use 5a s = 0.002 as the standard deviation. For ABM11 we compute symmetric 
errors with the formula 

Scr = ^Y^( a o -o-fc) 2 , (3.1) 

where the sum runs over all eigenvectors of the given set. This is the same prescription 
as in [IT]. Note that in the case of ABM 11 a s is one of the parameters of the fit. Thus, 
the PDF error computed above already includes the uncertainty due to the variation of 
the strong coupling constant. For the NNLL calculation we also have to determine the 
parameter (3 cut . The values are given in Table [5j 
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The results for the cross section at the Tevatron and at the LHC with ^/s = 7 TeV and 
y/s = 8 TeV are given in Table l6lR We find good agreement between the different PDF 
sets for the cross section at the Tevatron. At the LHC, we find again agreement between 
MSTW2008, NNPDF2.1, and CT10, but the value of ABM11 is significantly lower than 
the others. This reflects the smaller gluon density in the ABM11 fit that does not make use 
of the Tevatron jet data, which drive the high-x gluon distribution upwards. The difference 
in the top cross section is significantly larger than the error assigned to the individual PDF 
fits, which indicates that it may be possible to use the top-pair production cross section to 
constrain the gluon density. We pursue this further in Section [5j 

4 Heavy "top" quarks 

Heavy colour-triplet fermions appear in the extension of the standard model by a fourth 
generation, and are also often introduced to cancel the top-loop contributions to the Higgs 
mass in non-supersymmetric models that aim at addressing the hierarchy problem (such 
as Little Higgs models [3U[35]). In this context the extra fermions can be singlets under 
the SU(2)l gauge group or may transform under non-chiral representations, so they are 
less constrained than fourth generation quarks. Heavy colour-triplet fermions in other 
SU(2)l representations with exotic electric charges have also been introduced in models of 
electroweak symmetry breaking [36]. Within QCD, the production cross sections of these 
postulated particles only depend on the spin and colour charge and can be computed with 
Topixs. 

The LHC experiments have already performed direct searches for the production of 
heavy quarks in various decay channels, with the most stringent bounds mg > 611 GeV 
(for a heavy bottom quark B decaying exclusively as B — > tW, using 4.7 fb _1 of data [37]). 
rriT > 557 GeV (for a heavy top quark decaying exclusively as T — > bW, using 5.0 fb _1 
of data [38]), rriT > 475 GeV (for the decay T — > tZ allowed for SU(2) singlet T, using 
1.14 fb" 1 [31]). 

In Table [7J we present our predictions for the production cross section of heavy quarks Q 
at the LHC with 7 TeV centre-of-mass energy for the mass range trig = 350-700 GeV that 
is relevant for the analysis of 5 fb _1 of data collected in the 2011 run, and for masses up 
to 1 TeV at the LHC operating at y/s = 8 TeV. For each mass the value of the parameter 
Pcut appearing in the momentum-space resummation formalism with a running soft scale 
has been determined using the procedure introduced in p]. The resulting values range 
from f3 cut = 0.50-0.42 for m Q = 350-700 GeV and y/s = 7 TeV, and from f3 cut = 0.51-0.39 
for itiq = 350-1000 GeV at the LHC with ^/s = 8 TeV, and are implemented as default 
values in Topixs. We use the MSTW2008 PDFs and the running strong coupling with five 
active quark flavours and therefore neglect the impact of the top quark on the evolution of 
these quantities. The partonic cross sections are also computed with five massless flavours, 

8 The small difference in the MSTW2008 values with respect to Section [2] is due to the different value 
of a s (Mz). 
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Table 7: Total cross sections for the production of heavy quarks Q beyond the standard model 
in pb at the LHC (y/s = 7TeV, upper table and yfs = 8TeV, lower table). See caption of 
TableQ]for the definition of the various errors. Heavy "top" mass mg in GeV. 

LHC, x/s = 7TeV 



m Q NLQ NNLO app NNLL 2 



350 


3 


i i +0.32+0.17 
ii -0. 40-0.15 








3 


1 7+0.13+0.05+0.18 
■ l ' -0.16-0.05-0.16 








3 


91 +0.12+0.05+0.18 
• zi -0.11-0. 05-0.16 








400 


1 


07+O.I3+O.O8 
°' -0.17-0.07 








1. 


OQ+0.05+0.02+0.08 
OJ -0.07-0.02-0.07 








1 


a 1 +0.06+0.02+0.08 
• 4:i -0.05-0.02-0.07 








450 


6. 


11 +0.62+0.40 
•^-0.81-0.32 


X 


10 


-1 


6. 


[-£+0.23+0.09+0.40 
du -0.30-0.09-0.34 


X 


10 


-1 


6 


££+0.27+0.09+0.41 
uu -0.24-0.09-0.34 


X 


10 


-1 


500 


3. 


9n +0.30+0.21 
zu -0.40-0.16 


X 


10 


-1 


3. 


9^+0.11+0.04+0.21 
zu -0.14-0.04-0.17 


X 


10 


-1 


3 


O9+0.14+0.04+0.21 
• oz -0. 11-0. 04-0.17 


X 


10 


-1 


550 


1. 


fi7 +0.15+0.11 
u ' -0.21-0.09 


X 


10 


-1 


1. 


70+0.05+0.02+0.11 
■ ' u -0.07-0.02-0.09 


X 


10 


-1 


1 


70+O.O7+O.O2+O.II 
• ' °-0.06-0.02-0.09 


X 


10 


-1 


600 


8. 


Q7+0.82+0.63 
J ' -1.12-0.49 


X 


10 


-2 


9. 


1/1+0.26+0.09+0.63 
■ -^-0.38-0.09-0.48 


X 


10 


-2 


9 


01 +0.40+0.09+0.64 
• Oi -0.32-0.09-0.50 


X 


10 


-2 


650 


4. 


Q £+0.45+0.36 
JU -0.62-0.28 


X 


10 


-2 


5. 


n«+0.14+0.04+0.37 
uu -0.20-0.04-0.27 


X 


10 


-2 


5 


1 7+0.23+0.04+0.38 
1 ' -0.18-0.04-0.28 


X 


10 


-2 


700 


2. 


01 +0.25+0.22 
Oi -0.35-0.17 


X 


10- 


-2 


2. 


Q7+0.08+0.02+0.22 
■°' -0.11-0.02-0.16 


X 


10 


-2 


2 


qa +0.13+0.02+0.23 
• J ^-0.10-0.02-0.16 


X 


10 


-2 



LHC , y/s = 8 TeV 

m Q NLQ NNLO a pp NNLL^" 



350 


4 


nq+0.50+0.26 
JO -0. 62-0.23 








5 


n o+0.22+0.09+0.27 
uo -0. 26-0.09-0.24 








5 


n Q+0.19+0.09+0.28 
■ uy -0.18-0.09-0.24 








400 


2 


9O+0.22+0.12 
• Z6 -0.28-0.11 








2. 


97+0.09+0.04+0.13 
z '-0.ll-0.04-0.il 








2 


on+0-09+0.04+0.13 
' 3U -0.08-0.04-0.11 








450 


1. 


no+0.10+0.06 
UO -0.13-0.05 








1. 


-in+0.04+0.02+0.06 
iu -0.05-0.02-0.05 








1 


1 9+0.04+0.02+0.06 
iz -0.04-0.02-0.06 








500 


5. 


t-tr+0.51+0.33 
JJ -0.68-0.27 


X 


10 


-1 


5. 


fit -+0.19+0.07+0.34 
UJ -0.25-0.07-0.28 


X 10 


-1 


5 


7O+0.23+0.07+0.34 
• ' °-0.20-0.07-0.29 


X 


10 


-1 


550 


2. 


Q7+0.27+0.19 
J ' -0.36-0.15 


X 


10 


-1 


3. 


n o+0.10+0.04+0.19 
UO -0.13-0.04-0.15 


X 


10 


-1 


3 


n7 +0. 13+0.04+0.19 
• u ' -0.11-0.04-0.16 


X 


10 


-1 


600 


1. 


rk+0.15+0.11 
ud -0.20-0.08 


X 


10 


-1 


1. 


^Q+0.05+0.02+0.11 
UO -0.07-0.02-0.09 


X 


10 


-1 


1 


71 +0.07+0.02+0.11 
• ' i -0.06-0.02-0.09 


X 


10 


-1 


650 


9. 


a a +084 +0.63 
^-1.15-0.50 


X 


10 


-2 


9. 


£9+0.27+0.09+0.64 
uz -0.39-0.09-0.49 


x 10 


-2 


9 


70+0. 41+0.09+0.65 
• ' °-0.33-0.09-0.51 


X 


10 


-2 


700 


5. 


c 1+0.49+0.39 
d ^-0.68-0.30 


X 


10 


-2 


5. 


£1+0.16+0.05+0.39 
u ^-0.22-0.05-0.29 


X 


10 


-2 


5 


71+0.25+0.05+0.40 
'^-0.19-0.05-0.30 


X 


10 


-2 


750 


3. 


01 +0.29+0.24 
oi -0.41 -0.19 


X 


10 


-2 


3. 


qo+0.09+0.03+0.25 
oo -0.13-0.03-0.18 


X 


10 


-2 


3 


1 1+0.15+0.03+0.25 
•^-0.12-0.03-0.18 


X 


10 


-2 


800 


2. 


n9 +0.18+0.15 
uz -0.25-0.12 


X 


10 


-2 


2. 


n£+0.05+0.02+0.16 
uu -0.08-0.02-0.11 


X 


10 


-2 


2 


1 n +0. 09+0.02 +0.16 
iu -0.07-0.02-0.11 


X 


10 


-2 


850 


1. 


QK+0.11+0.10 
zd -0.15-0.08 


X 
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-3 
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10 


-3 
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O9+0.15+0.02+0.32 
• oz -0. 11-0. 02-0.21 


X 
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which corresponds to replacing the top quark t by Q (rather than adding Q to the known 
quarks). The same approximation was also used in [4]. 

For the mass range m Q = 350-700 GeV at the LHC with y/s = 7 TeV, the NNLL 
corrections grow from about 3% to above 4% relative to the NLO cross section evaluated 
with the appropriate NLO PDFs, while for tuq = 1 TeV at the LHC with y/s = 8 TeV 
the NNLL corrections become larger than 5%. In this comparison, the higher-order cor- 
rections are partially compensated by the change from NLO to NNLO PDFs, with the 
associated change in the value of the strong coupling. Normalizing the higher-order results 
instead to the NLO cross section evaluated with MSTW2008NNLO PDFs, the NNLL cor- 
rections range from 10% for itlq = 350 GeV to 14% for itlq = 1000 GeV at the LHC with 
y/s = 8 TeV, while the approximate NNLO cross section amounts to a 9-11% correction. 
Therefore higher-order effects beyond NNLO are small, but reach a magnitude of the order 
of the residual uncertainty of the approximate NNLO result for larger masses. 

The total theoretical uncertainty defined below ( 12.41) is largely independent of the mass 
and centre-of-mass energy. At yfs = 8 TeV, in the mass range vtlq = 350-1000 GeV, it 
ranges from to + ™% at NL0 > from -f% to -f% at NNLO, and from to + _fi at 
NNLL. It is counterintuitive that for larger masses the error estimate of the approximate 
NNLO prediction becomes smaller than that for NNLL. However, the size of our uncertainty 
estimate for large masses is similar to that found in |4j, so this issue is common to the 
Mellin-space and momentum space results. The relative PDF uncertainty is similar at 
NLO, NNLO a pp and NNLL and grows from tfi to for m Q = 350-1000 GeV at 
yfs = 8 TeV. Performing the computation with the NNPDF2.1 PDF set we find that the 
results for NNLL 2 (NNLO app ) are smaller by 4-6% (3%) if the same coupling constant is 
used for both PDF sets. Therefore both sets agree within the estimated PDF uncertainty. 

Since the contribution of the threshold region to the total cross section is more impor- 
tant for larger masses, we expect that the ambiguities in the implementation of threshold 
resummation are reduced and the predictions of different formalisms show better agree- 
ment than for lighter quarks. We therefore compare our default implementation, NNLL2, 
to an implementation NNLLfj x of resummation in momentum space with a fixed soft scale, 
determined using the method proposed in [TjJ], that was denoted as "Method 1" in [TJ. This 
implementation is available in Topixs with the option MUSRUN=0, and the default values 
of the fixed soft scale are implemented for MSTW2008 with m Q = 173.3-2000 GeV. By 
NNLLtv we denote the results obtained using the program top ++1.2 [5] that implements 
resummation using the traditional Mellin- moment method [4] . In Figure [1] the results of 
the three implementations are shown for the LHC with 8 TeV centre-of-mass energy for the 
mass range itlq = 200-1600 GeV. The larger masses are not relevant for phenomenology 
with 20 fb _1 luminosity, but are included here in order to compare the three implementa- 
tions in a regime where the threshold approximation is expected to be more accurate. We 
plot the quantity -K'nnll = cr tf NLL / cr tf LO ) where for comparison the NLO cross section in 
the denominator is evaluated with the NNLO PDFs, in contrast to our default treatment. 
It is seen that the three implementations are in good agreement once the higher-order 
corrections become larger than 10%, with larger differences for smaller masses. The dif- 
ference of the Mellin-space implementation to our default is always of a similar magnitude 
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Figure 1: Comparison of the ratio Knnll = CT jf NL V a tf L0 f° r heavy quark production at the 
LHC with y/s = 8 TeV for our default implementation NNLL 2 (black, solid) to our fixed- 
scale implementation NNLLfj x (blue, dashed) and the results NNLL^ (red, dot-dashed) from 
the Mellin-space approach [4][5]- The NNLL results are normalized by the NLO cross section 
computed with NNLO PDFs. 



or smaller than that of the two momentum-space implementations^ so we conclude that 
the momentum-space and Mellin-space approaches agree within the ambiguities inherent 
in the threshold approximation. 



5 Gluon distribution function 

The precise knowledge of parton distribution functions (PDFs) and of their uncertainties is 
a crucial requirement for the LHC physics program, particularly when the error on PDFs 
represents the largest single theoretical uncertainty, as is the case in top-pair production. 
PDFs are extracted by several collaborations [T0 l[TTll3Tll32ll40ll4T] from a large number 
of measurements collected at different experiments, the early fixed-target deep-inelastic 
scattering (DIS) and Drell-Yan (DY) experiments and the modern HERA and Tevatron 
colliders. In the meantime the LHC, thanks to the accumulated high statistics and the 
well-controlled systematics, is providing the PDF fitting collaborations with a large amount 
of precise data. Soon the LHC measurements are going to provide essential constraints 
on most PDF combinations. The NNPDF collaboration has already published a PDF 
analysis, NNPDF2.2, which includes early electro- weak LHC data from the 2010 runs [To] . 

9 The results from Topixs include bound-state effects and the resummation of Coulomb corrections 
beyond NNLO not included in top++. However, the bound-state contribution to Xnnll is of the order of 
0.3-0.4 percent for the mass range considered in Figure [TJ so the conclusion of the comparison to top++ 
is not affected. For top-pair production the contribution of the Coulomb corrections beyond NNLO above 
threshold was found to be even smaller [T], which we also expect to be the case for larger masses. 
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A recent study of the impact of LHC W lepton charge-asymmetry data on a MSTW2008 
Monte Carlo parton set was presented in [42j. 




In this section we investigate the possibility to constrain the gluon distribution in the 
proton through the top-pair production cross section by combining the most accurate 
available theoretical calculation presented in this paper with LHC data. This is interesting 
in several respects. First, while the gluon distribution function at medium-small x is 
rather well constrained by the precise HERA measurements of the scaling of the structure 
functions, at medium-large x the gluon is currently constrained primarily from the jet 
inclusive measurements performed at the Tevatron. A more precise determination of the 
gluon at medium-large x, which is the region probed by the ti measurements at the LHC, 
would lead to an improved theoretical precision of the Higgs cross section in gluon fusion. 
Moreover the top-pair production cross section is correlated to the strong coupling constant 
and therefore potentially reduces the uncertainty in the determination of a s from PDF 
global analyses. 

In this section, we quantify the impact of the LHC measurements of the total ti cross 
section on the gluon uncertainty by using the re- weighting technique [T^j[T5] . This technique 
is valid for Monte Carlo sets of PDFs, like those provided by the NNPDF collaboration. 
These sets determine a Monte Carlo representation of the probability density in the space of 
PDFs. This feature allows to include new experimental data by using Bayes' theorem, i.e. 
by re-weighting an existing PDF set (prior probability) without having to perform a new 
fit of the entire data. The effect of a new independent dataset is quantified by computing 
for each replica k (k = l,...,N Tep ) of the starting PDF set a weight w^, which assesses 
the probability that the PDF replica agrees with the new data. The re-weighted ensemble 
then forms a representation of the probability distribution of PDFs conditional on both 
the old and the new data. The weights are computed straightforwardly by evaluating the 
xt of the new data to each of the replicas, according to Eq. (9) of [TJ]. The distribution 
can then be unweighted and a new PDF set, which now includes both the data fitted in 
the initial PDF set and the additional data, is extracted. 

Before presenting the results of this analysis, it is important to specify the three in- 
gredients that are employed in order to determine the \\ °f ^ ne LHC top-pair production 
data to each replica and therefore the weights w^. the experimental data, the theoretical 
predictions and the prior PDF set. 

The (single) experimental data added to the PDF analysis is obtained from the com- 
bination of the ATLAS and CMS ti cross section measurements presented in [26 jH3l - H5] . 
While the proper combination of these results can only be performed by the experimental 
collaborations, we simply combine them here in a single data point by adding in quadrature 
statistical and systematic errors of the ATLAS and CMS combined results, and by adding 
the luminosity errors as a 100% correlated uncertainty between the two experiments, which 
results in 



a ex P . = 173 23 ± g 55 ^ stat + sygt ^ ± 7 0Q ( lumi _) pb = 173 _23 ± 9.59 p b. (5.1) 
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The theoretical predictions a^' k are obtained for each PDF replica k from (I2.2p l^l Since 
the PDF errors do not account for theoretical uncertainties, the latter are neither included 



in the computation of the weights Wk- As starting PDF sets, we consider two cases. In the 
first part of the analysis, the prior probability is given by the NNPDF2.1 NNLO parton 
set. The latter is a NNLO determination of parton distributions from a global set of hard- 
scattering data, including DIS, Drell-Yan and jet data. Each set is determined for a given 
value of a s (Mz), which is treated as an external parameter. Several sets are available with 
aP s = 0.114, 0.115, 0.124. Since the ti cross section is strongly correlated to a s , which in 
turn is strongly correlated to the gluon density, it is important to take into account the 
a s uncertainty in the prior PDF set. The N 3 ag replicas determining the starting parton set 
(such that Y^j=i N 3 as = 350) are taken from various a{ sets, as it was done in Section [3] 
to determine the PDF+a s uncertainty, according to a Gaussian distribution centered on 
a s (M z ) = 0.118 with a standard deviation of 5a s = 0.0015. 

The results of the re- weighting are displayed in Figure [2J On the left-hand side the im- 
pact of including the LHC measurement on the gluon density, with /if = m t = 173.3 GeV, 
is shown. As expected, there is an effect in the large- X region. The uncertainty in the gluon 
density given by the width of the bands can be assessed directly from the bottom-left plot, 
where the two lines refer to before and after including the ti cross section in the PDF 
analysis. Since the ti cross section predicted with the NNPDF2.1 PDF is close to (15.11) . 
see Table [6l there is only a small effect on the central gluon distribution. However, the 
uncertainty is reduced by about 20% for x > 0.15. The same results are shown on the 
right-hand side of Figure [2] but for the gluon-gluon luminosity, defined by 



where s is the LHC centre-of-mass energy squared and /i, following Ref. [16] . is set to 
fif = \frs. A similar impact as the one observed on the gluon density is observed on 
the gg luminosity in the expected region yfr > 2m t /y/s ~ 0.05 and partially also in the 
smaller yfr region due to the momentum conservation sum rules. The predicted ti cross 
section evaluated with the NNPDF2.1 set changes from crjr = 166.1 ± 7.0 (PDF) pb to 
ajj = 168.4 ± 5.4 (PDF) after re- weighting; the central value gets closer to the experimental 
data, and the PDF uncertainty is reduced by about 25%. On top of that, we observe that 
the x 2 P er degree of freedom to the Tevatron jet data is reduced by the inclusion of the ti 
measurement. If the same analysis were performed by using as a prior set the NNPDF2.1 
set with fixed a s , the impact would be rather smaller, as the reduction of the uncertainty 
is not only due to the gluon, but also to a s {Mz\ whose value computed over the replicas 
increases from 0.1180 ± 0.0015 to 0.1184 ± 0.0013 after including the ti data in the fit. 

It is interesting to observe that the effect is already visible with the data collected in 
the 2011 7 TeV LHC run. We can then expect an even more significant effect, when in 

ln The present analysis has been performed before [9] appeared and hence does not include the full qq 
NNLO partonic cross section. However, as discussed above, this has a negligible effect on the central value 
of the gluon-dominated ti cross section at LHC. Since the present re- weighting analysis does not include 
theoretical errors, only the central value matters. 




(5.2) 
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Figure 2: Effect of including the total ti cross-section measurement into the NNPDF2.1 global 
analysis for varying a s . Top: Gluon distribution (left) and gg luminosity (right) before (solid, 
solid green in color) and after (shaded, shaded blue in color) including the ti measurement in 
the NNPDF fit. Both curves are normalized to the central value of the NNPDF2.1 gluon before 
re-weighting. Bottom: Half-width of the error bands in the upper plots before (solid, solid 
green in color) and after (dashed, dashed blue in color) the inclusion of the ti measurement in 
the fit. 

the near future the experimental precision will reach the present theoretical accuracy. We 
therefore show in Figure [3] what happens if we sup po se that the measured ti cross section 
had an error as large as the theoretical uncertainty] 11 ! while the central value is kept fixed, 

alt = 173.231*3 pb, (5.3) 

and observe a further reduction by 5% (to a total of 25%) of the uncertainty in the gluon 
distribution to the NNPDF2.1 gluon in the large- x region. This means that, on the one 
hand, the limiting factor will soon be the theoretical uncertainty, which will be reduced 
further only once the full NNLO gg partonic cross section is known. On the other hand, 
it is remarkable to see a significant effect from a single additional measurement that com- 
petes with thousands of other measurements already included in PDF global analyses, in 

11 Here we adopt the error from the MSTW2008 PDF set in Table M that leads to a slightly smaller 
theoretical uncertainty. 
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Figure 3: Effect of including the total tt cross section measurement with reduced uncertainty 
into the NNPDF2.1 global analysis for varying a s . Gluon distribution (left) and (right) half- 
width of the gluon error band in the left plot before (solid light gray, solid green in color) and 
after including the tt measurement in the fit. In this case the dashed curve (dashed, dashed 
blue in color) refers to the inclusion of the real measurement (15. ID , the solid dark gray curve 
(solid, solid red in color) refers to the inclusion of the measurement with reduced uncertainty 
Q. 

particular with the approximately 200 inclusive jet data points collected at Tevatron Run 
II that are relevant in the large-x region. A larger effect should therefore be expected, once 
tt rapidity and invariant mass distributions will be included in the PDF fits. 

In order to disentangle the effect of the jet data from the one of the top-pair production 
data, we now start from the NNPDF2.1_DIS+DY set as prior probability, which has the 
same features as the NNPDF2.1 global set, but does not contain the jet data. This implies 
that the PDF uncertainty of the gluon at medium-large x is larger. We can therefore assess 
whether the top cross section pushes the gluon into the same direction as the Tevatron jet 
data and compare the constraining power of top-pair production relative to the inclusive 
jet measurements. Furthermore, the starting PDF set is supposedly closer to the ABM11 
NNLO PDF set [11], which does not include the Tevatron data, allowing us to draw some 
conclusion on the differences observed in Section |3]between theoretical predictions obtained 
with ABM11 and the other PDF sets. 

The theoretical prediction for the total top-pair production cross section obtained with 
the NNPDF2.1_DIS+DY set is a§- = 161.1 ± 9.9 (PDF) pb, 3% smaller than the one 
obtained with the NNPDF2.1 global analysis, but compatible within the PDF uncertainty. 
The x 2 P er degree of freedom of the tt data point to oQ' is 1.6. After including the tt 
measurement in the fit the predicted cross section moves to o^' = 167.5 ± 6.8 (PDF) pb 
and the % 2 decreases to 0.4. In FigureHJwe show the change of the gluon PDF and of the gg 
luminosity due to the new fit. The shift of the central value of the gluon is now significant, 
and it goes in the same direction as the shift when the jet data is included, bringing the 
black central line closer to one, which is the reference gluon of the global NNPDF2.1 fit. 
This can be expected since, as we observed earlier, the \ 2 P er degree of freedom to the 
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Figure 4: Effect of including the total tt cross section measurement (but excluding jet data) 
on the NNPDF2.1_DIS+DY PDF sets for varying a s . Top: Gluon distribution (left) and gg 
luminosity (right) before (solid, solid green in color) and after (shaded, shaded blue in color) the 
inclusion of the tt measurement in the NNPDF fit. Both curves are normalized to the central 
value of the NNPDF2.1 gluon. Bottom: Half-width of the error bands before (solid, solid green 
in color) and after (dashed, dashed blue in color) including the tt measurement in the fit. 

Tevatron jet data in the NNPDF2.1 global fit slightly improved with the inclusion of the 
tt measurement in the fit. A similar observation holds for the gg luminosity. The error 
reduction is significant and amounts to about 35% both for the gluon at large- x and the 
gluon-gluon luminosity at medium-large yjr. If we compare the size of the blue error band 
in Figure H] to the one of the green error band in Figure El we see that the error reduction 
due to the inclusion of the jet data is certainly larger, by a factor 4, due to the larger 
number of data included. However, it is striking that the single tt data point already has 
a non-negligible effect with respect to the Tevatron jet data. 

In Section[3]we observed that the MSTW2008, NNPDF2.1, and CT10 PDF sets predict 
comparable values for the top-pair production cross sections at the LHC, while the ABM 11 
prediction is significantly lower than the others. On the left-hand side of Figure [5] we com- 
pare the ABM11 gluon luminosity to the NNPDF2.1_DIS+DY and NNPDF2.1 ones, using 
the same a s (Mz) = 0.118 for consistency. In the large-x region the ABM11 gluon lumi- 
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GG luminosity normalized to NNPDF21 

NNPDF2.1,u 5 = 0.118± 0.0015 
| | NNPDF2.1(DIS+DY),ct s = 0.118±0.0015 

CD 





Figure 5: Comparison of gg luminosity obtained in the NNPDF2.1, NNPDF2.1_DIS+DY 
and ABM11 parton fits before (left) and after (right) the tt measurement is included in the 
NNPDF2.1 and NNPDF2.1_DIS+DY fits. All curves are normalized to the central value of the 
NNPDF2.1 gluon and computed for a s (M z ) = 0.118. 



nosity is about 2a smaller than NNPDF2.1. It is compatible with NNPDF2.1_DIS+DY, 
mainly because the uncertainty of the latter is much larger, which in turn is mostly due to 
differences in the gluon parameterization, and partially also due to a different treatment 
of the higher-twist effects in DIS. 

On the right-hand side of Figure [5] we show the same results but after including the 
top-pair production cross section in the two NNPDF2.1 fits. We already noted that this 
increases the large- x gluon and therefore the medium-large yjr gluon luminosity. A similar 
effect is expected on the ABM 11 gluon and should bring it closer to the other global 
determinations. But, due to the small uncertainty of the ABM11 gluon, the effect would 
be much more significant, and the inclusion of the it data might even require a more flexible 
large- X gluon parameterization. 

A caveat to the analysis in the present section must be mentioned, namely that it was 
performed for a fixed top-quark mass value of m t = 173.3 GeV. A change of ±1 GeV in the 
value of m t , which corresponds to the present top- mass uncertainty, changes the predicted 
inclusive top production cross section by roughly ±3%, which decreases the significance of 
this measurement for the NNPDF2.1 fit. On the other hand, the difference between the 
predicted it cross section obtained with the NNPDF2.1 set, = 166.1 ± 7.3 (PDF) pb, 
and with the ABM11 set, of- = 148.2 ± 5.9 (PDF)pb is about 10%, see Table El Hence, 
our conclusion that the present ti cross section measurement should already have an impact 
on the ABM11 fit remains valid after including the top-quark mass uncertainty. 



6 Conclusion 

In this paper we extended our analysis of the resummed total top-quark pair production 
cross section [TJ in several directions. We updated our results for Tevatron, and for LHC 
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with 7 and 8 TeV centre-of-mass energy to account for the full qq NNLO partonic cross 
section [9]; studied the dependence on different parton-distribution inputs and the impact 
of the LHC cross section measurement on the gluon distribution; and presented predictions 
for heavy "top" quarks. All results and many variations can be obtained with the publicly 
available program Topixs, which includes soft and Coulomb resummation. The main 
results of the present analysis can be summarized as follows: 

• Including the full qq NNLO partonic cross section removes the slight tension that 
existed in the prediction of the tt cross section at the Tevatron from different groups 
in favour of the somewhat larger values and reduces the uncertainty by almost a 
factor of two to about 3%, see also [9]. Our NNLO+NNLL Tevatron tt production 
cross section result, assuming m t = 173.3 GeV and a s (M z ) = 0.1171, is 

a& = 7.15 toil (theory) Igjjjj (PDF+a.) pb (MSTW2008NNLO). (6.1) 

The cross section increases to 7.26 pb, when the "world average" value a s (Mz) = 
0.1180 is used. A top mass m t = 171.4^" 7 GeV is obtained from Tevatron data. 

• Our prediction for the tt production cross section at the LHC with \fs = 8 TeV, 
assuming m t = 173.3 GeV and a s (M z ) = 0.1171, is 

<Ttt = 231.8 (theory) Ijj;? (PDF+a s ) pb (MSTW2008NNLO). (6.2) 

The cross section increases to 236.5 pb, when the "world average" value a s {Mz) = 
0.1180 is used. 

• The 2011 LHC data from the \fs = 7 TeV run decrease the uncertainty of the gluon 
distribution at x > 0.15 by about 20%, using NNPDF2.1 as the reference gluon. This 
effect is quite remarkable, given that only a single data point is added to the fit. The 
precision of the data and theoretical prediction will soon raise the issue of defining a 
procedure for including consistently theoretical uncertainties into the PDF fits and 
their error estimates. 

• PDF fits including the jet data from the Tevatron agree very well in their predicted 
top cross section at the LHC. The PDF uncertainty is now the dominant theoretical 
uncertainty. The ABM11 fit predicts a 10% lower cross section due to its smaller 
large-x gluon distribution. Our study leads us to expect that the tt cross section 
should already have a significant impact on this fit, even when all theoretical uncer- 
tainties, including the top-quark mass input are accounted for. 
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Appendix: Functionality of Topixs 



In this section we describe the main functionality of the program Topixs. A detailed 
manual is available together with the program package at the URL 

http : / /users . ph . turn . de/t31sof tware/topixs/ 

Topixs uses hplog [UJ for the numerical evaluation of harmonic polylogarithms and 
quadpack [18] for numerical integrations. The PDF sets are included via the lhapdf [19] 
library. The determination of the auxiliary parameter (3 cut and the fixed soft scale, which 
are used in the computation of the resummed cross section, uses routines from the GNU 
Scientific Library [50] . 

The program package consists of several FORTRAN and CH — h programs which implement 
the calculation of the total cross section as well as that of the parameter /? cut and the 
fixed soft scale. The user interface consists of a bash shell script, called topixs, which 
automatically compiles and executes these programs. The user can change the program 
settings by editing a single text file. Examples of such a configuration file are included in 
the program package. 

In our approach the cross section depends either on the auxiliary parameter /3 cut or 
the fixed soft scale. While programs for their computation are provided, we have also 
implemented default values for several PDF sets as well as the values for heavy quarks 
discussed in Section HJ Therefore, it is usually not necessary to compute them, in which 
case the GNU Scientific Library is not required either. 

Several approximations can be chosen for the cross-section computation. A fixed-order 
calculation is possible to NLO or NNLO( app ) accuracy. In the latter case, the user can 
choose whether or not to include the full result for the qq channel by changing the value 
of the variable QQNNLOEXACT in the configuration file. The resummed cross section can be 
evaluated in all of the different approximations defined in [I], i.e. (N)NLLi and (N)NLL2 
with fixed soft scale, and (N)NLL2 with running soft scale. The option BOUNDS allows the 
user to decide whether or not to include the contribution from bound states, or even to 
compute only this contribution. 

Our default choice is to use a running soft scale (MUSRUN=l) and include the bound- 
state contribution (B0UNDS=1). The full result for the qq channel is included in the NNLO 
computation and the matching of NNLL 2 (QQNNL0EXACT=1). These settings will be used 
unless the variables are explicitly reset in the configuration file. 

The main parameters of the calculation are the top-quark mass, the collider type 
(proton-proton or proton-antiproton) and energy, and the PDF set to be used. All of 
them can be chosen by the user. The default settings are m t = 173.3 GeV (MT0P=173 . 3), 
LHC with y/s = 8 TeV (C0LLIDER=1 and SQRTS=8000), and MSTW2008 (PDFSET=MSTW08). 
There are four additional predefined PDF sets: ABM11, NNPDF2.1, JR09 HQ], and CT10. 
For these sets the computation of the PDF+a s error is completely automated and it is 
possible to vary the value of a s (Mz) (except for JR09 where no fits with different values of 
a s (M z ) are provided). Other sets can be used by specifying the grid file name. It is also 
possible to choose a particular member of a PDF set to be used in the calculation. This 
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can be used to vary a s with a user specified grid file or for a re-weighting computation as 
in Section O 

The theory error is determined by the variation of the various scales (hard, soft, factor- 
ization, resummation, and Coulomb) and other parameters like /3 cut (cf. [lj for a detailed 
explanation). The central values of these scales can be chosen by the user. The variation 
is done automatically by evaluating the cross section at several equidistant points between 
twice and half the central value. The number of points can also be chosen by the user. 
Alternatively, the user can vary the scales by hand by performing the calculation without 
error computation (N0ERR=1) for different values of the scales. This also makes it possible 
to evaluate the cross section for specifically chosen points. 

The computation time strongly depends on the choice of the PDF set, with ABM11 
being the fastest and NNPDF2.1 the slowest (this was also noted in [H]). On an AMD 
Phenom 9850 processor with 2.5 GHz, the computation of a single cross section point 
(without any error estimates, using the default settings) takes approximately 15 seconds 
for ABM11, 75 seconds for MSTW2008, and 6 minutes for NNPDF2.1. Using the default 
configuration file to perform the computation of the NNLO app and NNLL2 approximations 
with MSTW2008, including the full theory and PDF+a s error, requires the evaluation of 
several hundreds of cross section points and takes about ten hours. 
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